NASA— CR-1 99085 


/-//Urf l— 


6/o? 3 



DESIGN OF ETO PROPULSION TURBINE USING CFD ANALYSES 


F.J. de Jong, Y-T. Chan, and H.J. Gibeling 


ORiviifivl. 

COLOR ILLUSTS 


» tiH 5 , 


Scientific Research Associates, Inc. 
Glastonbury, CT 06033 


April 1995 


Final Report R95-9084-F 


Prepared Under 
Contract NAS8-38865 


for 


National Aeronautics and Space Administration 
George C. Marshall Space Flight Center 


(NASA-CR-199085 ) DESIGN OF ETO 
PROPULSION TURBINE USING CFO 
ANALYSES Final Report (Scientific 
Research Associates) 89 p 


N95-32915 
Unc 1 as 


G3/37 0061028 



TABLE OF CONTENTS 


1. INTRODUCTION 1 

2. ANALYSIS 2 

2. 1 Navier-Stokes Equations 2 

2.2 Numerical Solution Procedure 3 

2.3 Grid Generation 4 

2.4 Boundary Conditions 5 

2.5 Initial Conditions 5 

3. RESULTS 6 

3.1 Baseline GGOT 6 

3.2 GGOT with Mini-Shroud 7 

3.3 Effect of Tip Clearance 7 

4. CONCLUSIONS 9 

REFERENCES 10 

FIGURES 12 

APPENDIX: Baseline Design of the Gas Generator Oxidizer Turbine (GGOT) and Performance 
Predictions for the Associated Oxidizer Technology Turbine Rig (OTTR) 



LIST OF FIGURES 


1 . Oxidizer Turbine Flow Path 

2. Baseline Turbine Blade Geometry 

3. Hub Section Grid 

4. Tip Section Grid (Baseline Geometry) 

5. Pressure Contours and Particle Traces oh the Baseline GGOT Blade at the Design 
Clearance 

6. Turbine Blade Geometry with Mini-Shroud 

7. Tip Section Grid (Geometry with Mini-Shroud) 

8. Pressure Contours and Particle Traces on the GGOT Blade with Mini-Shroud at the Design 
Clearance 

9. Blade Lift Coefficient as a Function of Radius 



1. INTRODUCTION 


State-of-the-art in pump design for space shuttle, space transport, or general ETO 
propulsion systems currently is a combination of experience, simple analyses with empiricism to 
estimate overall performance, and input from a database generated by experiments. This aspect of 
the design process will remain largely unchanged in the near future, due to the fact that current 
CFD viscous flow codes are "analysis" codes, rather than "inverse design" codes; i.e., they analyze 
the flow for a specified geometry and inflow conditions, rather than determine the geometry 
required to provide a desired flow field. With this limitation, CFD can be best utilized in modem 
pump design by first producing a "baseline" design produced by current design practice and 
utilizing state-of-the-art CFD codes to change design details so as to evolve the base design to an 
improved, advanced, and hopefully near-optimum design with improved performance. This process 
would utilize the insight which the computations provide into the flow field structure, to refine the 
baseline design or suggest new geometric configurations to achieve desired performance. 

In the work discussed here, a three-dimensional Navier-Stokes code was used for analysis 
of the STME fuel turbine, one of the tasks of the NASA/MSFC Turbine Technology Team (Refs. 
1-3). The STME baseline oxidizer turbine flow field was simulated and, based upon the simulation 
described here and simulations performed by other members of the NASA/MSFC Turbine 
Technology Team, an “advanced concept” design was developed by Pratt & Whitney, which, in 
turn, was analyzed. In addition, simulations were performed at different turbine blade tip 
clearances. 

The present report describes the CFD code used for these simulations, as well as results 
obtained. A report that details the design of the turbine by Pratt & Whitney has been included as an 
Appendix. 
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2. ANALYSIS 


2.1 Navier-Stokes Equations 

Solution of the flow field was obtained from a solution of the Reynolds-averaged, 
compressible Navier-Stokes equations. The governing equations were expressed in a rotating 
cylindrical coordinate system fixed to the turbine axis. In this coordinate system the mass, 
momentum, and energy conservation equations are: 


^ + v.( /J u)=0 

(1) 

^ • (pUU) +2pa)xU + pojxcox r = - V/>+ V* r 

(2) 

+ v*(pm) = Dp V«q + O 
at ’ Dt M 

(3) 


where U is the velocity vector in the rotating frame of reference, eo is the rotation vector, and r is a 
vector from the axis of rotation to the point under consideration. 

The stress tensor (molecular and turbulent) r is given by 

T ij = 2 ~ 2 A“<# V • U Sy (4) 


where the rate of the strain e tJ is given by 


«„=I/2 


dx j 


(5) 


and where the effective viscosity /^is the sum of the molecular and turbulent viscosities 


Mef =V + Vt ( 6 ) 

Here the turbulent viscosity pr is obtained from the turbulence model. is the viscous dissipation 
per unit volume, which can be expressed as 
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( 7 ) 


while the heat flux vector q is given by 

q = -(K+K r )Vr (8) 

Here tc and k t are molecular and turbulent thermal conductivities, respectively. In the present 
analysis, at and k t are obtained assuming constant molecular and turbulent Prandtl number Pr and 
Pr r , i e.. 


.t*p 

(9a) 

Pr 

_Mt c p 

Pr r 

(9b) 


A simple mixing-length type eddy viscosity model was used in the turbine computations, in 
which a mixing length distribution is specified normalized by a local freestream mixing length and 
modified to account for near-wall damping. The local freestream mixing length is proportional to a 
local shear layer thickness, which can be computed from the solution, or, as was done in the present 
calculations, can be specified. All boundary layers were assumed to be turbulent. 

2.2 Numerical Solution Procedure 


A general non-orthogonal coordinate transformation to a body-fitted grid is used to handle 
complex geometries in the solution procedure. The governing equations are solved by a Linearized 
Block Implicit (LBI) scheme (Refs. 4-5). 

The method can be outlined as follows: the governing equations are replaced by an implicit 
time difference approximation, optionally a backward difference or Crank-Nicholson scheme (a 
backward time-difference scheme was used in the present application). Terms involving 
nonlinearities at the implicit time level are linearized by Taylor series expansion about the solution 
at the known time level, and spatial difference approximations are introduced. The result is a 
system of multidimensional coupled (but linear) difference equations for the dependent variables at 
the unknown or implicit time level. To solve these difference equations, the Douglas-Gunn 
procedure for generating alternating-direction implicit (ADI) splitting schemes is introduced in its 
natural extension to systems of partial differential equations. This ADI splitting technique leads to 
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systems of coupled linear difference equations having narrow block-banded matrix structures which 
can be solved efficiently by standard block-elimination methods. Details are given in Refs. 4-5. 

In the present application, three-point central differences are used in the transformed 
coordinate system, and artificial dissipation terms of the form 




( 10 ) 


are added to the governing equations for each coordinate direction j. The variable <f> corresponds to 
the velocity component U t for the x r direction momentum equation, the density p for the continuity 
equation, and the enthalpy h for the energy equation. The coefficient (/W/ is obtained from the 
relation 

pU j AXj<(l/a d )^p+(p art ) j ^ ( 11 ) 


where Axj is the grid spacing at the point in question, while p corresponds to the effective viscosity 
Peff for the momentm equation, PerfPr for the energy equation, and is zero for the continuity 
equation. The artificial dissipation coefficient <jj lies between 0 (no dissipation) and 0.5 (full 
artificial dissipation). 

Use of artificial dissipation tends to enhance the stability and convergence properties of the 
numerical solution procedure, but it also tends to reduce the accuracy of the solution, in particular 
on coarse grids. Therefore, the turbine calculations were performed with an initial value of the 
artificial dissipation coefficient oj = 0.5 (corresponding to full artificial dissipation). After the flow 
field had been established, this coefficient was reduced to a* = 0 . 2 . This reduction of artificial 
dissipation did not affect the flow field qualitatively, but it did decrease the pressure loss and 
increase the efficiency. 

2.3 Grid Generation 


An important component in a three-dimensional Navier-Stokes simulation is grid 
generation. Grid generation for the present study was accomplished as follows. First, the EAGLE 
code (Ref. 6) was used to generate a two-dimensional grid for a number of blade sections. Next, 
these grids were “stacked” to form a three-dimensional grid. The grid points on hub-to-casing grid 
lines were redistributed to cluster points near the hub and near the tip by using Oh’s method (Ref. 
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7), and the diffusing section downstream of the blade was constructed by rescaling the local radii. 
In the tip clearance region, EAGLE was used to generate a two-dimensional grid “inside” the blade 
section at he tip (and matching the “exterior” grid). This grid was then used at all radial locations 
from the blade tip to the end wall. The resulting grid generation procedure is efficient, because it 
uses EAGLE only to generate two-dimensional grids, and uses an algebraic procedure to construct 
the three-dimensional grid from these two-dimensional grids. 

2.4 Boundary Conditions 

The computational domain chosen for the turbine calculations consisted of one passage 
between two blades, appropriately extended upstream and downstream of these blades. On this 
domain, the physical boundary conditions used were as follows: 

(i) No-slip and adiabatic wall conditions were specified on all solid surfaces. 

(ii) Total pressure, total temperature, and flow angles were specified at the inflow boundary. 

(iii) Static pressure was specified at the outflow boundary. 

(iv) Periodicity conditions were applied in the circumferential direction in the sections upstream 
and downstream of the blades. 

These boundary conditions were augmented by the appropriate numerical boundary conditions, viz. 
zero pressure gradient on stationary solid surfaces, zero reduced pressure gradient on rotating solid 
surfaces, extrapolation of pressure at the inflow boundary, and extrapolation of velocities and 
temperature at the outflow boundary. 

2.5 Initial Conditions 


In the turbine calculations, a steady-state solution was sought in the rotating frame of 
reference. Therefore, the initial conditions applied to these calculations serve as an initial guess, 
and do not affect the converged steady-state solution (although they will, in general, affect the 
convergence history). In the present calculations, an initial guess was obtained by calculating a 
two-dimensional solution for the mid-span geometry, and by using this solution for all blade 
sections. No attempt was made to include the tip clearance regions and hub or casing boundary 
layer profiles in the initial guess. 
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3. RESULTS 


The work presented here was generated as part of SRA's effort under the NASA/MSFC 
Turbine Technology Team. Results of previous work, the analysis of Pratt & Whitney’s Generic 
Gas Generator Turbine (GGGT), can be found in Refs. 8-9. The design of the present Gas 
Generator Oxidizer Turbine (GGOT), Pratt & Whitney’s effort under a subcontract from SRA, has 
been detailed in the Appendix (see also Ref. 10.). Most of the results of the present work, the 
analysis of the GGOT, have been presented at NASA MSFC meetings (see also Refs. 1 1 and 12.) 
Details have not been repeated here. Instead, a description is given of the different cases run, with 
a selection of representative results. 

3.1 Baseline GGOT 


The oxidizer turbine flow path and the turbine blade geometry are shown in Figs. 1-2. 
Additional information can be found in the Appendix. The present calculations simulate the rotor 
blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration 
terms included in the momentum equations (cf. Section 2.1). The upstream computational 
boundary is located about one axial chord from the leading edge. The boundary conditions at this 
location have been determined by Pratt & Whitney, using an Euler analysis without the vanes to 
obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage 
analysis including the vanes. Inflow boundary layer profiles were then constructed assuming the 
skin friction coefficient at both the hub and the casing. The downstream computational boundary is 
located about one axial chord from the blade trailing edge, and the circumferentially-averaged static 
pressure at this location was also obtained from the Pratt & Whitney Euler analysis. 

The simulation was run with a grid containing about 216,000 grid points: 90 grid points in 
the streamwise direction, 60 grid points in the circumferential direction, and 40 grid points in the 
hub-to-casing direction (with 14 grid points in the tip clearance region). Figures 3 and 4 show the 
grid at the hub section and the tip section, respectively. The tip section grid contains 37x31 grid 
points in the “interior” of the blade section. Figure 5 shows the calculated pressure contours on the 
blade and particle traces in the tip clearance region. The vortical flow behavior of the fluid that 
passes through the clearance region and exits at the downstream edge of the gap is clearly visible. 
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3.2 GGOT with Mini-Shroud 


Results obtained for the 3-D baseline GGOT geometry and the full-scale design Reynolds 
number show a region of high loss near the casing, attributed to the formation of the tip vortex (cf. 
Fig. 5). In an effort to reduce the clearance flow losses, the mini-shroud concept was proposed by 
the Pratt & Whitney design team: a “lip” or “mini-shroud” is added to the blade tip by replacing 
the blade pressure side at the tip by a straight line. The “thickness” of the mini-shroud (in the radial 
direction) has been kept fairly small (0. 1”). The airfoil shape below the mini-shroud is unchanged. 
Figures 6 and 7 show the geometry of the blade with the mini-shroud and the corresponding tip 
section grid used in the present calculations. To allow use of the grid generation procedure 
described in Section 2.3, the mini-shroud has been faired smoothly into the part of the blade below 
the mini-shroud. Results of the calculations indicate that, at the design tip clearance, the mini- 
shroud does not significantly affect the flow in the passage region. Although the tip region flow is 
affected, the results are qualitatively the same as those obtained without the mini-shroud, as can be 
seen by comparing Fig. 8 with Fig. 5. 

3.3 Effect of Tip Clearance 

To further investigate the effect of the mini-shroud on the tip clearance flow, calculations 
were carried out for both geometries (i.e., with and without the mini-shroud) at two additional tip 
clearances: 0.03” (1.5% span) and 0.06” (3% span), corresponding to 2x and 4x the design tip 
clearance, respectively. All calculations show flow patterns like those in Figs. 5 and 8 (see also 
Ref. 12). The results of the calculations have been summarized in Fig. 9, which shows the 
distribution of the blade lift coefficient as a function of radius for all six cases. For the baseline 
geometry, the unloading of the blade tip with increasing tip clearance is clearly visible. For the 
mini-shroud geometry, this effect is less pronounced. At the largest clearance, however, the hub 
loading decreases significantly for the mini-shroud geometry. But, since the calculations were 
performed with less grid resolution in the hub region than in the tip region, the prediction of the 
hub flow may not be quantitatively correct. 

To quantify the differences between the six cases, efficiencies were calculated. The “ideal 
work” was computed from the exit-to-inlet total pressure ratio (using isentropic relations); the 
“actual work” was computed from the decrease of tangential momentum from inlet to exit. In 
these efficiency calculations, the “exit” was taken about one chord downstream of the blade. Its 
precise location did not significantly affect the results. The “actual work” was also computed from 
the torque on the blade (without taking into account the effect of shear stress); the difference with 
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the work calculated from the tangential momentum was less than 1% in all cases. Table 1 lists the 
efficiencies as differences from the efficiency of the baseline case (the baseline geometry with the 
design tip clearance of 0.015”). Clearly, the efficiency decreases as the tip clearance increases. The 
mini-shroud seems to improve the efficiency a little bit at the design clearance and twice the design 
clearance, but this “improvement” is well within the accuracy of the calculations. At the largest tip 
clearance, the efficiency of the mini-shroud geometry has dropped off significantly, possibly due to 
hub-effects (cf. Figure 9). 

The conclusion that can be drawn from the calculations is that the effect of the mini-shroud 
on the tip clearance flow and the associated losses is not as large as was originally anticipated. 


Tip Clearance 

Baseline 

Geometry with 


Geometry 

Mini-shroud 

0.015” 

0 

+0.1% 

0.03” 

-1.2% 

-0.9% 

0.06” 

-2.4% 

-3.3% 


Table 1. Efficiencies Relative to the Baseline Efficiency 



4. CONCLUSIONS 


Under the present effort, a three-dimensional Navier-Stokes code was used for various 
turbine flow field calculations. The emphasis was on the flow in the tip clearance region and the 
associated losses. Both a baseline turbine blade geometry and an “advanced-concept” blade 
geometry (with a “mini-shroud”) were analyzed at different tip clearances. The results show that a 
current state-of-the-art CFD code can be used as an analysis tool in turbine design and the 
development of advanced hardware concepts. 
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FIGURE 1 - OXIDIZER TURBINE FLOW PATH 




FIGURE 2 - BASELINE TURBINE BLADE GEOMETRY 



FIGURE 3 - HUB SECTION GRID 
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FIGURE 4 - TIP SECTION GRID (BASELINE GEOMETRY) 



FIGURE 5 - PRESSURE CONTOURS AND PARTICLE TRACES ON THE 
BASELINE GGOT BLADE AT THE DESIGN CLEARANCE 




FIGURE 6 - TURBINE BLADE GEOMETRY WITH MINI-SHROUD 



FIGURE 7 - TIP SECTION GRID (GEOMETRY WITH MINI-SHROUD) 



FIGURE 8 - PRESSURE CONTOURS AND PARTICLE TRACES ON THE 

GGOT BLADE WITH MINI-SHROUD AT THE DESIGN CLEARANCE 



FIGURE 9 - BLADE LIFT COEFFICIENT AS A FUNCTION OF RADIUS 
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Introduction 


The work described here is being sponsored by the NASA Marshall Space Flight Cen- 
ter (MSFC), under contract No. NAS8-38865. This report is intended to summarize 
work completed under a Pratt & Whitney (P&W) subcontract from Scientific Research 
Associates (SRA) No. S900084-91-0001, Task I. Task I is as follows: 

The subcontractor shall provide the aerodynamic design for a baseline turbine with 
operating characteristics generally consistent with the NLS oxidizer turbopump, and 
shall provide follow-on support for the mechanical design of the baseline turbine rig. 

The objectives of this effort are to ensure that the aerodynamic design intent of the 
baseline turbine is achieved in the hardware, and to provide supporting aerodynamic 
analyses for the mechanical design effort. These analyses will include calculation of 
aero loads on vanes and blades, prediction of turbine output torque characteristics, and 
prediction of flowpath gas state conditions (pressure, temperature, and flowrate) for all 
turbine operating conditions of interest. Design consultation and review services will 
be provided to ensure the mechanical design of the flowpath meets aerodynamic de- 
sign requirements. A report detailing predicted turbine rig performance, torque, and 
flowpath conditions will be delivered. 

Turbine Technology Team 

The aerodynamic design of a technology demonstrator gas generator turbine, hereafter 
referred to as the Gas Generator Oxidizer Turbine (GGOT), was performed within the 
Consortium for Computational Fluid Dynamics Application in Propulsion Technology, by 
the Turbine Technology Team (Figure 1 ). Pratt & Whitney’s responsibility was to 
generate a preliminary turbine design for team member analyses, and to incorporate 
the results of these analyses into a baseline turbine design, which is the subject of this 
report (Figure 2). Consortium analyses are continuing in the "advanced concepts" 
phase. Results from these analyses along with the results from testing of the baseline 
turbine will be incorporated into an advanced concept turbine design. 

Design Table 

The design of the GGOT is consistent with the requirements for the U.S. Governments 
National Launch System (NLS) Oxidizer Turbine Design Table Revision 8A (Figure 3). 

In order to maximize specific impulse, a gas generator cycle with minimized turbine 
flowrate was chosen, which leads to high values of turbine specific work. The require- 
ments for high specific work and minimum weight led to the design of a highly loaded 
single stage oxidizer turbine which utilizes inlet and exit volutes to provide optimum 
performance. 

Aerodynamic Design Process 

The aerodynamic design of the baseline GGOT was initiated with Pratt & Whitney's 
generation of a preliminary turbine flowpath elevation along with stator and rotor airfoil 



contours. The process used to produce this preliminary configuration started with 
turbine meanline flowpath parametric analyses to optimize overall geometry such as 
blade height, flowpath diameter, and average blading reaction (Figure 4). Next, 
streamline analyses were conducted to optimize the radial distribution of flow proper- 
ties. The results of the streamline studies were used to create first pass 2D airfoil 
sections at a number of spanwise locations for both stator and rotor. These sections 
were radially faired to build 3D models which were then analyzed in multi-stage 3D 
Euler and Navier-Stokes flow solvers. Numerous iterations involving airfoil contour 
refinement were performed to arrive at the preliminary design. 

Once the preliminary design was completed, Turbine Technology Team members 
selected various areas for more detailed analyses. These included a 2D unsteady 
Navier-Stokes analysis of the turbine stage and 3D Navier-Stokes analyses of the vane 
alone, blade alone, turbine stage, and blade tip leakage details. The results from these 
studies were then incorporated into the design to define the baseline turbine (Figure 5). 

Inlet Volute Design 

The aerodynamic design of the baseline inlet volute manifold was accomplished using 
a Computer Aided Design (CAD) program to generate trial contours, which were then 
analyzed with a 3D Navier-Stokes flow solver. The design features axircular cross 
section, with an approximately linear decrease in thru-flow area in the circumferential 
direction (Figures 6A-6C). A decreasing area volute was selected over a geometrically 
simpler constant area toroid for several reasons: predicted reduction in the inlet radial 
side load by 80%, reduction in weight, and the creation of a highly swirled and uniform 
exit flow which enables a reduction in stator count. 3D Navier-Stokes pressure distribu- 
tions indicate that the inlet volute should experience a minimal transverse static pres- 
sure gradient at the aerodynamic design point (Figure 7). 

Turbine Stage Design 

The baseline oxidizer turbine is a single stage configuration which rotates in a clock- 
wise direction when viewed from the aft looking upstream (Figure 8). This highly 
loaded turbine uses a mixture of oxygen and hydrogen gases to deliver over 1900 
horsepower at a flow rate of 60.4 Ibm/s. The turbine stator airfoil contour design has 
been integrated with that of the inlet volute (Figure 9A-9G). This integration of flowpath 
components has significantly reduced stator turning, and the number of stators re- 
quired. 3D Euler analysis predicts that the stator will smoothly accelerate the flow from 
a Mach number of 0.35 to 1 .05 through a combination of minimal turning and endwall 
convergence at the aerodynamic design point (Figure 10A-10E). 

In an effort to improve performance while reducing cost, this turbine has been designed 
with reduced thru-flow velocities and increased blade reaction, relative to previous gas 
generator turbine experience. This causes the required blade turning angle to increase 
to an average of 157 degrees (Figure 1 1 A-1 1G). While this level of gas turning is 
beyond traditional design experience, 3D Euler analysis predicts that the rotor will 
smoothly accelerate the flow from a Mach number of 0.80 to 1.05 (Figure 12A-12E). 



Numerous independent CFD analyses performed by Turbine Technology Team mem- 
bers have consistently confirmed the prediction of smooth airfoil velocity distributions. 
This design concept is predicted to improve turbine stage efficiency by as much as 2%, 
while reducing the required blade count by 50% relative to traditional designs. 

Exit Volute Design 

The aerodynamic design of the baseline exit volute manifold was accomplished using 
the same methods and similar criteria as the inlet volute manifold. However, the chal- 
lenge of designing the exit volute is compounded by the higher Mach number and 
diffusing flowfield. To ensure turbine operation over the entire test matrix, two exit 
volutes were designed. One exit volute was designed with an oversized square cross- 
section to ensure that the flow did not choke at off-design operation (Figure 13A-13C). 
The other exit volute was designed with a round cross-section to smoothly collect the 
turbine exit flow at the aerodynamic design point and provide optimum efficiency (Fig- 
ure 14A-14C). The exit volute is positioned 139 degrees counter-clockwise relative to 
the inlet volute when viewed from the top of the engine configuration looking down 
(Figure 15). 

The turbine exhaust system consists of three regions: an annular diffuser located 
downstream of the rotor, a volute manifold with an increasing flow area, and a conical 
pipe diffuser for static pressure recovery (Figure 16). The annular diffuser reduces the 
Mach number by approximately 5% while acting as a short buffer zone between the 
turbine and exit volute. The volute flowpath has an approximately linear increase in 
area in the circumferential direction to minimize the static pressure gradient. A conical 
diffuser, located at the exit volute exit reduces the Mach number to a level required by 
the downstream ducting and engine nozzle. 

3D Navier-Stokes pressure distributions indicate that the square oversized exit volute 
will experience a substantial static pressure gradient (Figure 17A-17B). However, this 
design will allow the turbine to operate over the entire test matrix, dropping the flow 
Mach number from 0.80 to 0.60 at the aerodynamic design point. 3D Navier-Stokes 
pressure distributions indicate that the round, aerodynamically designed, exit volute will 
experience a minimal transverse static pressure gradient (Figure 18A-18B). 

Oxidizer Technology Turbine Rig 

The Oxidizer Technology Turbine Rig (OTTR, Figure 19) is a 50% scale model of the 
GGOT previously described. The OTTR, like the hot engine design, is an integrated 
turbine system consisting of an inlet vclute, a turbine stage, and an exhaust volute and 
diffuser. Aerodynamic performance evaluation of the OTTR will be performed in 
MSFC's Turbine Test Equipment (TTE) facility. 

Similarity requirements were used to define the "cold" turbine test article operating 
conditions (Figure 20). Testing at reduced inlet pressure and temperature facilitates 
the collection of quality data, but requires the scaling of other parameters to maintain 
flow similarity. Scaling was done as follows: 


Pressure Ratio 


( Ptin / Ptex ) 


Speed Parameter ( N D / V z R Ttin ) 

Flow Parameter ( 4 W V z R Ttin / n D 2 Pt, N ) 

Where, Ptin Inlet Total Pressure 
Ptex Exit Total Pressure 

N Shaft Speed 

D Rotor Mean Diameter 
z Compressibility 

R Gas Constant 

Ttin Inlet Total Temperature 
W Weight Flow 

Pressure ratio and speed parameter will be set when testing the OTTR, while flow 
parameter will be an output. Reynolds number, another similarity parameter will not be 
matched but evaluated as follows: 

Reynolds Number : (PsVX/zRTsp) 

Where, Ps Local Static Pressure 

V Local Velocity 

X Airfoil Axial Chord 

Ts Local Temperature 

p Fluid Viscosity 

Reynolds number is predicted to be a secondary influence on turbine performance 
(Figure 21 ). To evaluate this effect, the following turbulent boundary layer Reynolds 
number correlation can be used (Glassman): 

Correlation: ( 1 -Tl) GGOT 

d-TDon* 

Where, r| Turbine Efficiency 

Rn Average of Vane & Blade Exit 
Reynolds Numbers 

The value of the exponent was derived from Pratt & Whitney's extensive turbine test 
data base. 

OTTR Performance Predictions 

OTTR test objectives include measurement of turbine performance at the aerodynamic 




design point (ADP) and at off-design conditions (Figure 22). Nominal inlet conditions 
for testing of this turbine will be 100 psia and 560 R. Design point pressure ratio and 
speed at these inlet conditions will be 1.60 and 3710 rpm, respectively. The off-design 
envelope will include pressure ratios from 1 .20 to 1.80, and speeds from 1000 to 5000 

rpm. 

Turbine performance, such as efficiency and flow characteristics, will be mapped and 
compared to predictions (Figures 23 and 24). Other performance parameters such as 
flowpath pressures, temperatures, Mach numbers, and air angles, as well as turbine 
power output will also be assessed (Figure 25A-25E). 

The performance of the total turbine system will be evaluated and compared to design 
intent. Losses will be measured for the various components to determine where design 
goals are being accomplished and where improvements can be made (Figure 26A- 
26B). Also, detailed flow measurements including laser velocimitry will be used to pin- 
point areas of potential improvement. 

Summary 

The design of the Gas Generator Oxidizer Turbine (GGOT), and the associated test 
article, the Oxidizer Technology Turbine Rig (OTTR) are being done in support of 
turbine research and development within the Turbine Technology Team, which is part 
of the Consortium for Computational Fluid Dynamics Application in Propulsion Technol- 
ogy The baseline OTTR will undergo performance evaluation testing in MSFC's short 
duration experimental turbine test facility. The baseline test results will be integrated 
with analytical results of ongoing team studies to define an "advanced concept oxidizer 
turbine system having potential for further performance and cost benefits. This turbine 
design is the result of a successful team working relationship, in which organizations 
with various turbine backgrounds have worked together in a complementary and syner- 
gistic manner. Turbine geometric data sets (including the inlet volute, vane, blade, and 
exit volute) may be obtained via the Turbine Technology Team Coordinator, Mail Stop 
ED-32, National Aeronautics and Space Administration, George C. Marshall Space 
Flight Center, MSFC, AL 35812. 



FIGURE I 

OXIDIZER TURBINE BASELINE DESIGN 

Consortium for CFD Application in Propulsion Technology 
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FIGURE 2 . 

OXIDIZER TURBINE BASELINE DESIGN 

Turbine Technology Team - Oxidizer Turbine Baseline Design Completed 



Hot Fire Testing 















FIGURE 3 

OXIDIZER TURBINE BASELINE DESIGN 

Hot T/P Design Requirements (NLS Design Table Revision 8A) 
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FIGURE 4 

OXIDIZER TURBINE BASELINE DESIGN 

Aerodynamic Design Process 



Turbine Geometry 






FIGURE 5 

OXIDIZER TURBINE BASELINE DESIGN 

Full Scale Turbine Flowpath 



Axial Location - ( ) 


FIGURE £>A 

OXIDIZER TURBINE BASELINE DESIGN 

Inlet Volute - 50% Scale OTTR Flowpath 




Axial Location - ( ) 



FIGURE £>B 

OXIDIZER TURBINE BASELINE DESIGN 

Inlet Volute - Unigraphics Wire Frame Mesh 







FIGURE &G 

OXIDIZER TURBINE BASELINE DESIGN 

Inlet Volute - Geometry Definition 
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FIGURE 7 

OXIDIZER TURBINE BASELINE DESIGN 

Inlet Volute - Predicted Static Pressure Distribution at Vane Inlet Plane 



Circ. Position 


FIGURE 8 

OXIDIZER TURBINE BASELINE DESIGN 

Single Stage Oxidizer Turbine - 3D Computational Grid Viewed from Front 






FIGURE 9A 

OXIDIZER TURBINE BASELINE DESIGN 

Vane Elevation 



NOTES: 1 . Reference Drawing 80M52760 For Machining Details 

2. Number of Vanes is 20. 

3. Vane is Curved Line Faired Through Defining Sections. 




FIGURE 9B 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Stacked Airfoil Sections 





FIGURE 9C 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Section at Radius=4.4865 




FIGURE 9D 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Section at Radius=4.7489 





FIGURE 9E 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Section at Radius=5.0113 





FIGURE 9F 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Section at Radius=5.2736 




FIGURE 9G 

OXIDIZER TURBINE BASELINE DESIGN 







FIGURE IOA 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Pressure Distribution at 0% Span 




FIGURE 10B 

OXIDIZER TURBINE BASELINE DESIGN 



FIGURE IOC 

OXIDIZER TURBINE BASELINE DESIGN 





FIGURE IOD 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Pressure Distribution at 75% Span 
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FIGURE IOE 

OXIDIZER TURBINE BASELINE DESIGN 

Vane - Pressure Distribution at 100% Span 
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FIGURE II A 

OXIDIZER TURBINE BASELINE DESIGN 

Blade Elevation 



NOTES: 1 . Reference Drawing 80M52752 For Machining Details 

2. Number of Blades is 42 

3. Blade is Curve Line Faired Through Defining Sections 

4. Rotation is Clockwise when viewed from Aft looking Forward 




FIGURE IIB 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Stacked Airfoil Sections 





FIGURE IIC 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Section at Radius=4.4865 
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FIGURE l ID 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Section at Radius=4.7489 





FIGURE HE 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Section at Radius=5.0113 





FIGURE IIF 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Section at Radius-5.2736 





FIGURE IIG 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Section at Radius=5.5360 





FIGURE 12 A 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Pressure Distribution at 0% Span 
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FIGURE 12 B 

OXIDIZER TURBINE BASELINE DESIGN 

Blade - Pressure Distribution at 25% Soan 





FIGURE I2C 

OXIDIZER TURBINE BASELINE DESIGN 




FIGURE I2D 

OXIDIZER TURBINE BASELINE DESIGN 





FIGURE I2E 

OXIDIZER TURBINE BASELINE DESIGN 





FIGURE I3A 

OXIDIZER TURBINE BASELINE DESIGN 

Square Exit Volute - 50% Scale OTTR Flowpath 




FIGURE 136 

OXIDIZER TURBINE BASELINE DESIGN 

Square Exit Volute - Unigraphics Wire Frame Mesh 







FIGURE I3C 

OXIDIZER TURBINE BASELINE DESIGN 

Square Exit Volute - Geometry Definition 
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FIGURE 14 A 

OXIDIZER TURBINE BASELINE DESIGN 

Round Exit Volute - 50% Scale OTTR Flowpath 
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FIGURE 146 

OXIDIZER TURBINE BASELINE DESIGN 

Round Exit Volute - Unigraphics Wire Frame Mesh 
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FIGURE I4C 

OXIDIZER TURBINE BASELINE DESIGN 

Round Exit Volute - Geometry Definition 
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FIGURE 15 

OXIDIZER TURBINE BASELINE DESIGN 

Inlet Volute & Exit Volute Relative Position 
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FIGURE \(o 

OXIDIZER TURBINE BASELINE DESIGN 




FIGURE I7A 

OXIDIZER TURBINE BASELINE DESIGN 



FIGURE 176 

OXIDIZER TURBINE BASELINE DESIGN 

Square Exit Volute - Predicted Flow Patterns 
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FIGURE I8A 

OXIDIZER TURBINE BASELINE DESIGN 




FIGURE 18 B 
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FIGURE 19 

OXIDIZER TURBINE BASELINE DESIGN 

Oxidizer Technology Turbine Rig (OTTR) Flowpath (50% Scale of Hot T/P) 



Axial Location - ( ) 


FIGURE 20 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR Design Point and Operating Envelope 
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FIGURE 21 

OXIDIZER TURBINE BASELINE DESIGN 

Predicted Reynolds Number Effect on Efficiency 
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FIGURE 11 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR Test Matrix 



Speed (rpm x 1 000) 




FIGURE 23 

OXIDIZER TURBINE BASELINE DESIGN 

0777 ? - Predicted Efficiency 
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FIGURE 24 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Predicted Flow Parameter 
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FIGURE 25A 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Predicted Performance at Ptin=100 psia and Ttin=560 °R 
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FIGURE 25B 

OXIDIZER TURBINE BASELINE DESIGN 
















FIGURE 25C 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Predicted Performance at Ptin=100 psia and Ttin=56 
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FIGURE 25D 

OXIDIZER TURBINE BASELINE DESIGN 












FIGURE 25E 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Predicted Performance at Ptin=100 psia and Ttii 
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FIGURE 2(oh 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Component Performance with Square Exit Volute 


























FIGURE 2L 6 

OXIDIZER TURBINE BASELINE DESIGN 

OTTR - Component Performance with Round Exit Volute 
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